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The relativistic mean-field approach was implemented in the Lanzhou quantum molecular dynamics trans- 
port model (LQMD.RMF). Using the LQMD.RMF, the properties of collective flow and pion production were 
investigated systematically for nuclear reactions with various isospin asymmetries. The directed and elliptic 
flows of the LQMD.RMEF are able to describe the experimental data of STAR Collaboration. The directed flow 
difference between free neutrons and protons was associated with the stiffness of the symmetry energy, that is, 
a softer symmetry energy led to a larger flow difference. For various collision energies, the ratio between the 
a and r” yields increased with a decrease in the slope parameter of the symmetry energy. When the collision 
energy was 270 MeV/nucleon, the single ratio of the pion transverse momentum spectra also increased with 
decreasing slope parameter of the symmetry energy in both nearly symmetric and neutron-rich systems. How- 
ever, it is difficult to constrain the stiffness of the symmetry energy with the double ratio because of the lack of 


threshold energy correction on the pion production. 
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I. INTRODUCTION 


The equation of state (EOS) of nuclear matter, which orig- 
inates from the interactions between nuclear matter, plays an 
important role in nuclear physics and astrophysics. Heavy- 
ion collisions, the properties of nuclei, and neutron stars 
(NSs) have been widely studied to extract the nuclear EOS. 
Because nuclear many-body problems are highly nonlinear 
and the EOS is not a directly observable quantity in experi- 
ments, there are still some uncertainties in the EOS despite 
great efforts [1-6]. For instance, the EOS extracted from 
the GW170817 event has uncertainties at high nuclear den- 
sities [2]. Although the EOS can be extracted from the prop- 
erties of NSs, the internal composition of NSs is still poorly 
understood. The core of an NS may contain exotic mate- 
rials such as hyperons, kaons, pions, and deconfined quark 
matter [7—12]. Heavy-ion collisions in terrestrial laboratories 
provide a unique opportunity to study both the EOS and ex- 
otic materials. 

Collective flows of heavy-ion collisions were proposed in 
the 1970s and first detected in experiments at Bevalac [13- 
16]. Because collective flows are associated with nucleon- 
nucleon interactions, nucleon-nucleon scattering, etc., collec- 
tive flows have been used to extract the nuclear EOS [1]. Col- 
lective flows are also helpful for understanding the phase tran- 
sition between hadronic and quark matter. Generally, when a 
phase transition between hadronic and quark matter occurs, 
the collective flows of heavy-ion collisions indicate a soft 
EOS [17-21]. In addition, the ratios of the isospin parti- 
cles in heavy-ion collisions, such as 7~/a+, K°/K*, and 
=~ /*, are thought to be sensitive to the stiffness of the 
EOS [22-28]. The production of pions and kaons has been 
experimentally measured in 197 Au+!° Au collisions. The 
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K* production predicted by various transport models favors a 
soft EOS of isospin-symmetric nuclear matter at high baryon 
densities [29-33]. The m~ /r* ratio predicted by various 
transport models is still model-dependent [34-37]. Based 
on the FOPI data for the n~ /n* ratio [38], some results fa- 
vor a stiff symmetry energy (isospin asymmetric part of the 
EOS) [34, 35], whereas others imply a soft symmetry en- 
ergy [36, 37]. Recently, by analyzing the ratios of charged 
pions in 197Sn +124 Sn, 12Sn +14 Sn, and 19%Sn +"? Sn 
collisions [39], a slope of the symmetry energy ranging from 
42 to 117 MeV was suggested [40, 41]. The collective flows 
and ratios of charged pions are still worth studying to find the 
sources of the difference in various transport models and to 
extract information about the EOS from heavy-ion collisions. 

Quantum molecular dynamics (QMD) is a popular trans- 
port model that has been developed into many versions 
and used to successfully describe the properties of nu- 
clear matter, nuclei, mesons, and hyperons [33, 34, 42-55]. 
In high-energy heavy-ion collisions, the relativistic effects 
should be considered in QMD because they become signif- 
icant. The RQMD approach has been proposed for this pur- 
pose [42, 43]. Recently, relativistic mean field theory (RMF) 
was implemented in RQMD (RQMD.RMF) [44-46]. The 
RQMD.RMEF has been successfully applied to investigate the 
collective flows of hadrons [44—46]. In this study, we im- 
plemented RMF theory with isovector-vector and isovector- 
scalar fields in the Lanzhou quantum molecular dynamics 
model (LQMD.RMF). The channels for the generation and 
decay of resonances (A(1232), N*(1440), N*(1535), etc.), 
hyperons, and mesons were included in the previous LQMD 
model [33, 34, 47-50]. Using the LQMD.RMF, we explored 
the relationship between the symmetry energy and the prop- 
erties of the collective flow and pion production. 

The remainder of this paper is organized as follows. In 
Sect. II, we briefly introduce the formulas and approaches 
used in this study. The formulas include RMF theory, the 
dispersion relation, and the production of pions. The results 
and discussion are presented in Sect. III. Finally, a summary 
is presented in Sect. IV. 


I. FORMALISM 
A. Relativistic mean field approach 


The RMF interaction is achieved by exchanging mesons. 
Scalar and vector mesons provide medium-range attrac- 
tion and short-range repulsion between nucleons, respec- 
tively [56]. The nonlinear self-interaction of the ø meson is 
introduced to reduce the incompressibility to a reasonable do- 
main [57]. To investigate the properties of symmetry energy, 
we also consider the isovector-vector p [58] and isovector- 
scalar ô mesons [59]. The Lagrangian density is expressed 
as [59, 60] 
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where My = 938 MeV is the nucleon mass in free space, gi 
with i = g, w, p, ô is the coupling constant between nucleons, 
m; with i = 0,w,p,6 denotes the meson mass, gz and g3 are 
the coupling constants for the nonlinear self-interaction of the 
a meson, and Fy = pwy — O,w,, and Buy = Oyby — Osby 
are the strength tensors of the w and p mesons, respectively. 
The equations of motion for the nucleon and meson are ob- 
tained from the Euler-Lagrange equations and are written as: 
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where p and ps are the baryon and scalar densities, respec- 
tively, 3 = Pp — Pn is the difference between the proton and 
neutron densities, and ps3 = PSp — 08x is the difference be- 
tween the proton and neutron scalar densities. 

In the RMF approximation, the energy density is given by 
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where pr denotes the nucleon Fermi momentum, and Mž = 
My — goo F 9563 (— proton, + neutron) denotes the effec- 
tive nucleon mass. The isospin splitting of the effective nu- 
cleon mass My — M% still has a large uncertainty at this point. 


Analyses of nucleon-nucleus scattering data based on the op- 
tical model favor My — My > 0[61, 62]. Calculations based 
on Brueckner theory [63—65] and density-dependent relativis- 
tic Hartree-Fock [66] also indicate M% — M; > 0. How- 
ever, Mă — Mě < 0 is predicted by the transport model 
for heavy-ion collisions [67, 68] and nonlinear RMF mod- 
els [59, 60, 69]. In addition, both M% — M; < 0 and 
M, — M% > 0 can be predicted by the point-coupling RMF 
[69] and Skyrme and Gogny forces [70-76]. Because the La- 
grangian density in this study is the same as that in Ref. [59] 
and [60] except for the parameter settings, as shown in Fig. 1, 
the negative isospin splitting of the effective nucleon mass 
Mă — Mg < 0 is consistent with that in Ref. [59] and [60]. 
In the nonlinear RMF model, the isospin splitting of the ef- 
fective nucleon mass is primarily determined by the coupling 
strength of the 6 meson. The large coupling strength of the 
ô meson results in large isospin splitting of the effective nu- 
cleon mass. When the coupling strength of the 6 meson is 
zero (setl), there is no isospin splitting of the effective nu- 
cleon mass. 
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Fig. 1. (Color online) Difference between neutron and proton effec- 
tive masses as a function of the baryon density. 


Using the isospin asymmetry parameter a = (fn — 
Pp) /(Pn + Pp), the symmetry energy can be written as [59] 
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where f; = i, i = p, ô, and Ef = ypg + M*? and M* = 
Mn — 900 are the effective nucleon masses of the symmetric 
nuclear matter. The integral A(pp, M*) is defined as 
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Fig. 2. (Color online) Symmetry energy as a function of the baryon 
density. 


can be obtained as 
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In this study, we set the saturation density to pọ = 0.16fm~°, 
and the binding energy per particle of the symmetric nuclear 
matter was set to E/A — My = —16 MeV. For symmetric 
nuclear matter, we set setl, set2, and set3 models to be the 
same as the result of vanishing isospin asymmetry. As shown 
in Table 1 and Fig. 2, the symmetry energy of setl, set2, and 
set3 was set to be 31.6 MeV at the saturation density. Set1 
contained only p mesons; however, set2 and set3 contained 
both the p and 6 mesons. For setl, when the symmetry energy 
was set to be 31.6 MeV at the saturation density, the coupling 
parameter g, was fixed and the slope of symmetry was fixed 
at L = 85.3 MeV. For set2 and set3, when the symmetry 
energy was fixed at 31.6 MeV, the slope of the symmetry en- 
ergy was obtained by varying the coupling parameters g, and 
gs. Because the effective mass M* of the above models for 
symmetric nuclear matter was the same, the symmetry energy 
with both the p and ô mesons could not be softer than that of 
setl containing only p mesons. To broaden the range of the 
slope parameters, we set the slope parameter of set2 and set3 
to be 109.3 and 145.0 MeV by varying the coupling param- 
eters gp and gs, respectively. A broader range of slope pa- 
rameters would be helpful for understanding the relationship 
between the properties of symmetry energy and the observ- 
ables of heavy-ion collisions. The curvature of the symmetry 
energy Kym, which is a higher-order expansion coefficient 
of the symmetry energy compared to the slope parameter L, 
may also affect the properties of the symmetry energy and the 
observables of heavy-ion collisions at densities far beyond the 
saturation density. Kym of set1, set2, and set3 is obtained as 
—15, 141, and 391 MeV, respectively. 


B. Relativistic quantum molecular dynamics approach 


To investigate high-energy heavy-ion collisions, RQMD 
was proposed [42, 43]. Recently, RMF was implemented in 
RQMD [44-46]. In RQMD, for an N-body system, there 
are 4N position coordinates q} and 4N momentum coordi- 
nates p! (i = 1,..., N). However, the physical trajectories 
(q; and p;) are 6N for an N-body system. 2N constraints 
are required to reduce the number of dimensions from 8N to 
physical trajectories 6N [42-46, 77-79]. 


a (11) 


where 2N constraints satisfy the physical 6N phase space. 
The sign ~ indicates Dirac’s weak equality. The on-mass 
shell conditions can reduce the phase space from 8N to 7N 
dimensions. 
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where 2 = 1,..., N. The remaining N constraints are time fix- 
ation constraints. A simple choice of time fixation constraints 
that obey the worldline condition is written as [43, 77, 79, 80] 


Mj? = (p; - Vi)? - (My - 8)? =0, (12) 


Qin =4- (Gi — an) =0, (i = 1,...,N — 1), 

g2n =4-qn —T=O0, (13) 
where â = (1, 0) denotes a four-dimensional unit-vector [42— 
46, 77]. In a two-body center-of-mass system, â is defined as 
Pig / T with pi; = p} + p}. We observe that only the con- 
straint 2 = 2N depends on 7. With the above 2N constraints, 


the number of dimensions from 8N will reduce to 6N. These 
2N constraints are conserved over time. 
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Because only the constraint 1 = 2N depends on 7, A is writ- 
ten as[77] 
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with Cz = [¢;, d;]. Following previous studies, the Hamil- 
tonian of the N-body system was constructed as a linear com- 
bination of 2N — 1 constraints [77, 79, 80]: 
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Assuming [¢;,¢;] = 0, A; = 0 for N +1 <i < 2N [77]. 
The equations of motion are then obtained as 
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Table 1. Parameter sets for RMF. The saturation density po is set to 0.16 fm~*. The binding energy of the saturation density is E /A— My = 
—16 MeV. The isoscalar-vector w and isovector-vector p masses are fixed at their physical values, mu = 783 MeV and mp = 763 MeV, 
respectively. The remaining meson masses, Mo and ms, are set to be 550 and 500 MeV, respectively. 


Model go gu go(fm') g3 Ip g K(MeV) Esym(po)(MeV) L(po)(MeV) M*(po)/My Ksym (MeV) 
setl 8.145 7.570 31.820 28100 4.049  - 230 31.6 85.3 0.81 -15 
set? 8.145 7.570 31.820 28.100 8673 5347 230 31.6 109.3 0.81 141 
set3 8.145 7.570 31.820 28.100 11.768 7.752 230 31.6 145.0 0.81 391 


with the on-mass shell conditions (Eq. (12)) as inputs, the 
equations of motion can be obtained as 
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where z7” = p;"/p*° and M* = My — Si. The scalar 
potential S; and vector potential V;,, in RQMD are defined as 


Si = Goo + gstiði, 
Viu — Bi gu Win + Bitigpbi,u, (19) 


where t; = 1 for protons, and t; = —1 for neutrons, and B; 
is the baryon number of the ith particle. The meson field can 
be obtained from the RMF: 
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Substituting Eq. (20) into Eq. (19), the scalar potential S; 
and vector potential V;,, of the nucleons can be obtained. 
For other hadrons, such as A resonances, similar to other 
transport models [34, 85], the A optical potential is esti- 
mated using the nucleon potentials and square of the Cleb- 
sch—Gordan coefficient. In the RQMD approach, the scalar 
density, isovector-scalar density, baryon current, and isovetor 
baryon current are expressed as 
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Because the difference between the numerical results ob- 
tained using the effective mass M; and kinetic momentum 
ph * in the density and current and those obtained using a free 
mass M; = My = 938 MeV and canonical momentum p% in 
the density and current is small, a free mass M; = My = 938 
MeV and canonical momentum Di have been used in the 


above density and current [45]. The interaction density pij 
is given by a Gaussian 
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squared [77], and yi; is a Lorentz factor that ensures the 
correct normalization of the Gaussian [81] and is equal to 


(p? + p§)/ (pi + pj) in the two-body center-of-mass frame. 
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In this study, we set the Gaussian width to Lg = 2.0 fm?. 


C. Dispersion relation and production of pions 


The Hamiltonian of the mesons is defined as [48, 82-84] 


Nm 
Hy = X [VE + wi, ps]; (23) 


i=l 
where V£ is the Coulomb potential, which is expressed as 


DE eje; 
a (24) 
ja "9 
and Ny and Ng are the total numbers of mesons and baryons, 
including charged resonances, respectively. The pion poten- 
tial in the medium, which contains isoscalar and isovector 


contributions, is defined as 
w (Di, pi) = Wisoscalar (Pis pi) + C,T720(p/ po), (25) 


where a denotes the isospin asymmetry parameter, the co- 
efficient C, is 36 MeV, the isospin quantity 7 is 1, 0, and 
-1 for 7~, 7°, and at, respectively, and yy determines the 
isospin splitting of the pion potential and is set to two. In this 
study, the isoscalar part of pion potential Wisoscalar Was Chosen 
as the A-hole model. The pion potential, which contains a 
pion branch (smaller value) and A-hole (larger value) branch, 
is defined as 


n (Di, Pi Wr —tike (Di, Pi) + 


Wisoscalar (pi ’ pi) = 
Sa (Di, pi)WA-iike( Di, pi). (26) 


The probabilities of the pion and A-hole branches satisfy the 
following equation: 


The probability of both the pion and A-hole branches is de- 
fined as [84] 


1 


Sf) = r nra. 28 
where w denotes wr—like and WA-tike. The eigenvalues of 
Wr—like aNd wA-—iike are generated from the pion dispersion 


relation 
w? = P + mÈ +HI(w), (29) 


where II denotes the pion self-energy. Including the short- 
range A-hole interaction, the pion self-energy is defined as 


Lox 
where m, denotes pion mass. The Migdal parameter, g’, was 
set to 0.6. x is defined as 
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WA — 
where wa = Vmi + De — My and ma is the delta mass. 
In this study, the 7 NA coupling constant fa was 2, and the 
cutoff factor b was TMz. 

Both the Coulomb and pion potentials contribute to the de- 
cay of resonances and the reabsorption of pions. For instance, 
the energy balance of A° in the decay of resonances and the 
reabsorption of pions can be written as 


(30) 


X= 


Vm + oR = y MÈ + (Pr — Pr)? + wr (Prs P) + VE, 
(32) 


where pr and p, are the momenta of the resonances and pi- 
ons, respectively, and mp is the mass of the resonances. Be- 
cause A° is uncharged, the Coulomb potential exists only for 
charged particles after the decay of A°. 

The pion is generated from direct nucleon-nucleon colli- 
sion and decay of the resonances A(1232) and N*(1440). 
The reaction channels of the resonances and pions, which 
are the same as those in the LQMD model, are as fol- 
lows [33, 48, 88, 89]: 


NNONA, NN &NN’*, 
N* & Nr, 


For the production of the A(1232) and N* (1440) resonances 
in nucleon-nucleon scattering, the parameterized cross- 
sections calculated using the one-boson exchange model were 
employed [91]. 

The decay width of A(1232) and N*(1440), which orig- 
inates from the p-wave resonances, is momentum-dependent 
and is expressed as [91] 


a.\pl* 
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where |p] is the momentum of the created pion (in GeV/c). 
The parameters a1, a2, and a3 were taken as 22.48 (17.22) 
c/GeV, 39.69 (39.69)c?/GeV?, and 0.04(0.09) GeV?/c?, re- 
spectively, for A(N*). The bare decay width of A(N*) 
was given by T'o = 0.12(0.2)GeV. With the momentum- 
dependent decay width, the cross-section of pion-nucleon 
scattering has the Breit-Wigner form: 


7A 0.2512 
Onn(V/8) = Omax (2 (P) Q? (35) 
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where p and Pm are the three momenta of the pions at ener- 
gies of \/s and mọ, respectively. The maximum cross-section 
Omax Of A and N* resonances was obtained by fitting the 
total cross-sections of the experimental data in pion-nucleon 
scattering using the Breit-Wigner formula [92]. For instance, 
the maximum cross-section Omax Of A resonance was 200, 
133.33, and 66.7 mb for rtp > Att (mn > A), 
np + At (mn —> A®°) anda-~p— A} (ntn > At), 
respectively [89]. 

Note that the threshold effect was neglected in this 
study. The threshold effect mainly refers to the A pro- 
duction threshold energy and incident energy of two col- 
liding nucleons modified by the medium. The inci- 
dent and threshold energies in the medium are defined 


as \/Sin = VE + E? + Ex + 09)? — (X1 + Hg)? and 
sh = Vins 4+ E9 + mx + D2)? — (Š; + ¥E4)2, respec- 
tively [85-87]. Because Š; = 0 and p; ~ 0 for static nu- 
clear matter, the difference between the incident and thresh- 
old energies is \/Sin — \/Sm = Ef + ©} + ES + XQ — m3 — 
E9 — mă — X9 [86]. The difference between the incident and 
threshold energies, which is isospin dependent, may result in 


an enhancement in the nn — pA- channel and suppression 
of the pp + nAt* channel. 


Ill. RESULTS AND DISCUSSIONS 


The directed and elliptic flows were derived from the 
Fourier expansion of the azimuthal distribution: 


um) = Noll + 2V1 (y, pr)cos(¢) 
+2V2(y, pr)cos(2¢)], (36) 


where the azimuthal angle of the emitted particle ø was mea- 


sured from the reaction plane. pr = ,/p? + p3, denotes the 


transverse momentum, and the directed flow V; and elliptic 
flow V2 are expressed as follows: 


Vi = (cos(d)) = (2), 


Di 300 
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The directed flow provides information on the azimuthal 
anisotropy of the transverse emission. The elliptic flow de- 
scribes the competition between the in-plane (V2 > 0) and 
out-of-plane (V2 < 0) emissions. 

Firstly, the collective flows of LQMD.RMF in the !97 Au + 
197 Au collisions have been investigated at an incident energy 
of 2.92 GeV/nucleon (the corresponding nucleon-nucleon 
center-of-mass energy is vV Sny = 3 GeV ). The collec- 
tive flows of LQMD with Skyrme interaction and without 
the momentum-dependent interaction have also been inves- 
tigated. The Skyrme interaction of symmetric nuclear mat- 
ter is taken to be the same as that of SLy6, with an incom- 
pressibility of K = 230 MeV at pọ = 0.16fm~°[93, 94]. 
The symmetry energy of Skyrme interaction is defined as 


2 
Esm = 4 he (252)7/3 + S (p/po)™. When the Csym and 
Ys are set to be 38.6 MeV and 1.049, respectively, the sym- 
metry energy and the slope parameter of symmetry energy 
are 31.6 MeV and 85.3 MeV, respectively. The symmetry 
energy and the slope parameter of symmetry energy of the 
Skyrme interaction are as same as those of setl. As shown in 
Fig. 3, we have compared the collective flows of LQMD.RMF 
and LQMD with Skyrme interaction with recent experimental 
data from STAR Collaboration [90]. The collective follows 
of LQMD.RMF and LQMD with Skyrme interaction can de- 
scribe the STAR data at an impact parameter b =4 fm and 
b =7 fm, respectively. The directed flows of LQMD.RMF 
are almost consistent with the STAR data across the entire ra- 
pidity. However, The directed flows of LQMD with Skyrme 
interaction are weaker than the STAR data across the entire 
rapidity. This result may be due to the fact that the value of di- 
rected transverse momentum with the Lorentz effect is larger 
than that without the Lorentz effect [43]. The elliptic flows of 
both LQMD.RMF and LQMD with Skyrme interaction are 
consistent with the STAR data at rapidities smaller than 0.5. 
However, the elliptic flows of both LQMD.RMEF and LQMD 
with Skyrme interaction are weaker than the STAR data at 
large rapidity. At an incident energy of 2.92 GeV/nucleon, the 
LQMD.RMEF can better describe the experimental data com- 
pared to the LQMD with Skyrme interaction and without the 
momentum-dependent interaction. Based on the above analy- 
ses, the RMF models have been implemented into the LQMD 
model successfully. 


With this LQMD.RMF model, the }°°Sn + 1!2Sn and !8?Sn 
+ !24Sn collisions in this study were investigated at an inci- 
dent energy of 270 A MeV and an impact parameter b = 3 
fm. At an incident energy of 270 A MeV, the nuclear matter of 
the collision center can be compressed to densities approach- 
ing 2/9. In this dense region, collective flows, which reflect 
nucleon-nucleon interactions, can be used to extract the high- 
density behavior of the EOS [1, 44, 79, 88, 95]. The collective 
flows of free protons in the !°°Sn + 1!?Sn and !°?Sn + !74Sn 
collisions are shown in Fig. 4 and Fig. 5, respectively. It is 
reasonable that the maximum value of the directed flow V; 
was significantly larger than that of the elliptic flow V2. In 
the same reaction system, the difference in the directed flows 
with various slopes of symmetry energy (setl, set2, and set3) 
was small. The difference in the elliptic flows with various 
slopes of symmetry energy was also small. To determine the 
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Fig. 3. (Color online) Rapidity distribution of the collective flows 
of free protons in the 1°’ Au + +97 Au reaction at an incident energy 
of 2.92 GeV/nucleon. The open circles correspond to the STAR 
data [90]. 


relationship between the slope of the symmetry energy and 
the collective flow, we must process the collective flow data. 

The difference between neutron and proton directed flows 
emitted from heavy-ion collisions can be used to extract the 
density dependence of the symmetry energy [88, 95]. The dif- 
ference between the neutron and proton directed flows is de- 
fined as Vj, — Vip, as shown in Fig. 6. The trend and shape of 
the difference between the neutron and proton directed flows 
were similar to those of nonrelativistic LQMD [88]. For a 
given reaction system (nearly symmetric !°°Sn +11? Sn sys- 
tem or neutron-rich !°?Sn+!?4Sn system), the absolute value 
of the difference between the neutron and proton directed 
flows with a soft symmetry energy was higher than that with a 
stiff symmetry energy. Interestingly, the stiffness of the sym- 
metry energy can be reflected through the difference between 
the neutron and proton directed flows. The relationship be- 
tween the curvature of the symmetry energy Ksym and the col- 
lective flows was then investigated. The difference between 
Kym Of set] and Kym of set3 was 406 MeV, which was sig- 
nificantly larger than the 59.7 MeV difference between L of 
setl and L of set3. Although the curvature of the symmetry 
energy Kym also affected the difference between the neutron 
and proton directed flows, because the nuclear matter of the 
collision center could only approach 29 at an incident en- 
ergy of 270 A MeV, the effects of Kym were not significant 
compared to the effects of the slope parameter L. 

In addition to collective flows, the production of isospin 
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Fig. 4. (Color online) Rapidity distribution of the collective flows of 


free protons in the ‘°°Sn +"? Sn reaction at an incident energy of 
270 MeV/nucleon. 
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Fig. 5. (Color online) Rapidity distribution of the collective flows of 
free protons in the ®°Sn +1'** Sn reaction at an incident energy of 
270 MeV/nucleon. 
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Fig. 6. (Color online) Difference between neutron and proton di- 
rected flows in the }°8Sn +1!? Sn and +° Sn +14 Sn reactions at 
an incident energy of 270 MeV/nucleon. 
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Fig. 7. (Color online) Ratio between the 7 and m* yields as a 
function of the incident energy in the '°*Sn +14 Sn reaction. 


exotic particles such as hyperons, kaons, and pions can also 
be used to extract the symmetry energy [22-28]. Because the 
thresholds of hyperons and kaons were not reached at the in- 
cident energies in this study, the isospin exotic particles were 
pions. First, we calculated the ratio between the m~ and zt 
yields of the neutron-rich 1°?Sn +124 Sn system as a func- 
tion of the collision energy at the impact parameter b=3 fm 
and Om < 90°. Because setl had the softest symmetry en- 
ergy, it had the highest neutron density. Consequently, there 
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Fig. 8. (Color online) Transverse momentum spectra of pions as 
functions of transverse momentum at an incident energy of 270 
MeV/nucleon. The left two panels [(a) and (b)] are the results of 
the ?°°Sn +1? Sn reaction, and the right two panels [(c) and (d)] 
are the results of the 1??Sn +174 Sn reaction. The full circles corre- 
spond to the S7RIT data [40]. 


were more neutron-neutron scatterings in setl, resulting in 
the production of more A~ and 7~. As shown in Fig. 7, the 
m/n" ratio of setl was the highest, and the n~ /7* ratio 
of set2 was higher than that of set3. Specifically, at a colli- 
sion energy of 270 MeV/nucleon, the m~ /m* ratio without 
(with) the r potential changed from 2.71 (2.54) to 2.23 (2.06) 
when the slope parameter was varied from L =85.3 to 145.0 
Mey, that is, from setl to set3. In other words, the 17 / at 
ratio as a function of collision energy depends on the stiff- 
ness of the symmetry energy. This result was consistent with 
the predictions of most transport models [28, 39, 86]. When 
the r potential was considered, the interaction between 7 and 
the nucleon became attractive, resulting in a decrease in both 
n~ and 7+ via the absorbed channels 7N —> A(1232) and 
A(1232)N — NN. However, with the 7 potential, because 
there were more neutrons in the neutron-rich 82Sn +174 Sn 
system, 7~ was more easily absorbed than 7*. Consequently, 
the m~ /r* ratio without the m potential was higher than that 
with the 7 potential in the same RMF model. Moreover, it 
is worth mentioning that the threshold effect neglected in this 
study may cause the m~ /r* ratio to be reversed. In other 
words, with the threshold effect [85-87], 7 / a? of set3 may 
be the highest, and the m~ /r* ratio of set2 may be higher 
than that of setl. 

Next, the properties of m were predicted as a function of 
the transverse momentum. As shown in Fig. 8, the left and 
right panels are the transverse momentum spectra of pions 
for the nearly symmetric 1°8Sn +!!? Sn and neutron rich 
132Sn +124 Sn reactions at Am < 90°, respectively. In col- 
lisions between isotopes, 7+ is mainly generated from colli- 
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Fig. 9. (Color online) Single spectral ratios of pions as functions 
of transverse momentum for the !°°Sn +*!? Sn and '8?Sn +14 Sn 
reactions at an incident energy of 270 MeV/nucleon. The full circles 
correspond to the S7RIT data [40]. 


sions between protons and 77 is mainly generated from col- 
lisions between neutrons. Theoretically, a stiffer symmetry 
energy would have a stronger repulsive force to push out neu- 
trons and a stronger attractive force to squeeze protons, re- 
sulting in a decrease in the m~ yield and an increase in the 
mt yield, respectively. As shown in panels (b) and (d) of 
Fig. 8, the stiffer symmetry energy indeed led to larger trans- 
verse momentum spectra for 7+. However, the relationship 
between the transverse momentum spectra of m~ and the stiff- 
ness of the symmetry energy could not be directly explained. 
Compared with the stiffness of the symmetry energy, the 7 
potential had a more significant impact on the transverse mo- 
mentum spectrum of 7. For the neutron-rich !8?Sn +174 Sn 
system, the transverse momentum spectra of both 7+ and m7 
without the 7 potential were lower than those of the StRIT 
data at pr < 200 MeV. When the m potential was consid- 
ered, the transverse momentum spectra of both nt and m7 
increased at pr < 200 MeV but decreased at pr = 200 MeV. 
Consequently, the transverse momentum spectra of 7+ with 
the m potential were almost consistent with the S7RIT data 
[40]; however, the transverse momentum spectra of 7~ were 
still lower than the S7RIT data for the entire pr domain. The 
lower transverse momentum spectra of 7~ may be due to the 
absence of a threshold effect. The threshold effect, which was 
not considered in this study, may enhance the production of 
n~ [85-87]. 

Because a stiffer symmetry energy would have a stronger 
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Fig. 10. (Color online) Double ratio of pions as a function of trans- 
verse momentum at an incident energy of 270 MeV/nucleon. The 
full circles correspond to the S7RIT data [40]. 


repulsive force to push out neutrons and a stronger attractive 
force to squeeze protons, resulting in a decrease in the m7 
yield and an increase in the 7* yield, respectively, the single 
ratio SR(n—/n*) = [dM (x—)/dpr]/[dM(x*)/dpr] may 
depend on the stiffness of the symmetry energy and the reac- 
tion system. As shown in Fig. 9, for both the nearly symmet- 
ric system and neutron-rich system, a softer symmetry energy 
led to a larger single ratio. In addition, for the same stiffness 
of the symmetry energy, because the neutron-neutron scatter- 
ing of the neutron-rich '°?Sn+1*4Sn system was greater than 
that of the nearly symmetric }°°Sn +11? Sn system, the sin- 
gle ratio of the neutron-rich system was higher than that of 
the nearly symmetric system. It is worth mentioning that the 
single ratio of 1°°Sn +1! Sn was almost consistent with the 
experimental data. However, the single ratio of !?Sn+1!74Sn 
was lower than that of the experimental data for the entire pr 
domain. This was because the transverse momentum spectra 
a of 13?Sn +124 Sn were lower than the experimental data. 


The double ratio between the neutron rich system 
and the nearly symmetric system DR(a~/m*) = 
SR(a7 /7*)1394194/SR(a7 /t*) 1084112; which can 
cancel out most of the systematic errors caused by Coulomb 
and isoscalar interactions, was suggested to extract the 
properties of the symmetry energy [40]. However, because 
a lower symmetry energy will lead to a larger single ratio for 
both the nearly symmetric system and neutron-rich system, 
as shown in Fig. 10, it is still difficult to understand the 
dependence of the double ratio on the stiffness of symmetry 
energy. In addition, the double ratio without the 7 potential 
decreased slightly as a function of the transverse momentum, 
whereas it increased slightly as the transverse momentum 
increased when the m potential was considered. However, the 
increasing trend was still considerably weaker than that of 
the experimental results. The lower double ratio originated 
from the lower single ratio of the neutron-rich 1°?Sn +174 Sn 
system compared with the experimental data. The threshold 
effect may enhance the production of m~ [85-87] and the 


single ratio of the neutron-rich system, resulting in a larger 
double ratio. 


IV. CONCLUSION 


An RMF with various symmetry energies, namely setl, 
set2, and set3, was implemented in LQMD. The collective 
flows of the nearly symmetric '°°Sn +1! Sn and neutron- 
rich !32Sn +1%4 Sn systems were successfully generated from 
the LQMD.RMEF. It has been observed that the directed flow 
Vı was an order of magnitude larger than the elliptic flow 
V2. However, the difference between the directed flows V; 
with various slopes of symmetry energy was small. The dif- 
ference in the directed flows V2 with various slopes of sym- 
metry energy was also small. To explore the relationship be- 
tween the collective flow and the stiffness of the symmetry 
energy, we defined the difference between neutron and proton 
directed flows Vin — Vip. For a given nearly symmetric sys- 
tem or neutron-rich system, the absolute value of Vi, — Vip 
increased with decreasing slope of the symmetry energy. 

We also investigated the relationship between isospin ex- 
otic particles and the stiffness of the symmetry energy. The 
ratio between m~ yield and m* yield of the neutron-rich 
132Sn +124 Sn system as a function of the collision energy 
increased with a decrease in the slope parameter of the sym- 
metry energy. At an incident energy of 270 MeV/nucleon, the 
properties of 7 were predicted as a function of the transverse 
momentum. For the nearly symmetric !°°Sn +11? Sn system, 
the single ratio of the nearly symmetric system was consistent 
with the experimental data. However, for the neutron-rich 
132Sn +124 Sn system, the single ratio was lower than the 
experimental data, resulting in a double ratio lower than the 
experimental data. The 7 potential did not explain the lower 
transverse momentum spectra of 7~ in the neutron-rich sys- 
tem. Considering the 7 potential, the double ratio increased 
slightly with increasing transverse momentum. However, the 
increasing trend was still considerably weaker than that ob- 
served in the experimental results. The single ratio of the 
neutron-rich system and the double ratio may be lower than 
the experimental data because of the lack of a threshold ef- 
fect. The threshold effect, which can enhance the production 
of m~, could be a candidate for enhancing the single ratio of 
the neutron-rich system to a double ratio. Moreover, because 
a softer symmetry energy led to a larger single ratio for both 
nearly symmetric and neutron-rich systems, the dependence 
of the double ratio on the stiffness of the symmetry energy 
was not significant. The sensitivity of the double ratio to the 
stiffness of the symmetry energy may also be due to the lack 
of a threshold effect. When the threshold effect is considered, 
the production of 7~ in a neutron-rich system may be more 
significant than that in a nearly symmetric system. In the fu- 
ture, when collective flows, the single ratio of the neutron-rich 
system and the double ratio of the LQMD.RMF are consistent 
with the experimental data, Vi, — Vıp and the single ratio of 
the neutron-rich system m~ /r™, which are sensitive to the 
stiffness of the symmetry energy, may be used to extract the 
slope parameter of the symmetry energy. 


Appendix A: DETAILS of EQUATION OF MOTION 


For numerical calculations, the equation of motion (Eq. 
(18)) must be written in the computed form. With Eq. (19) 
and M¥ = M; — S; = My — Sj; as the inputs, the equation 
of motion (Eq. (18)) can be expanded as 
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In the two-body center-of-mass frame, p;; equals p;;. The 
squared distance q7 ;; is reduced to gz ;; = —7j; aa ii) 
, , z 
In the actual calculation, we replace p? with \/p? + M? to 
save calculation time [79]. Thus, the partial derivative of den- 
sity versus momentum and space can be written as 
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